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Abstract 

In this talk we examine the azimuthal angle decorrelation in dijet production at large 
rapidity intervals at the Tevatron as a possible signature of the BFKL evolution. 
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1 Introduction 



The state-of-the-art in jet physics at hadron coUiders is described by next-to-leading- 
order (NLO) QCD parton-level calculations [|], They appear to be in very good 
agreement with the one- and two-jet inclusive distributions obtained from the data of 
the CDF experiment at the Fermilab Tevatron Collider |Q, [Q. However, since data 
are being collected at the CDF and DO detectors at larger and larger rapidities, it 
may be possible to probe kinematic configurations where this fixed-order analysis is 
inadequate. This could occur when the cross section contains logarithms of large ratios 
of kinematic invariants. Typical invariants are the hadron-hadron center-of-mass energy 
y^, the parton-parton center-of-mass energy \/§ = ^/xaXbs, where xa and xb are the 
momentum fractions of the partons originating the hard scattering, and the momentum 
transfer Q, which is of the order of the transverse momentum of the jets produced in the 
hard scattering. If xa and xb are large, then the large logarithms, ln(s/Q^), factorize 
into the partonic subprocess cross section. These logarithms, which are of the size of the 
rapidity interval in the scattering process, can be resummed by using the techniques of 
Balitsky, Fadin, Kuraev, and Lipatov (BFKL)0. 

In analyzing dijet production experimentally so that it most closely resembles the 
configuration assumed in the BFKL theory, the jets are ordered first by their rapidity 
rather than by their energy]^. Thus, we look at all the jets in the event that are above a 
transverse momentum cutoff p±min, using some jet-definition algorithm, and rank them 
by their rapidity. We then tag the two jets with the largest and smallest rapidity and 
observe the distributions as a function of these two tagging jets. The cross section is 
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inclusive so that the distributions are affected by the hadronic activity in the rapidity 
interval y between the tagging jets, whether or not these hadrons pass the jet-selection 
criteria. We will refer to these hadrons in the rapidity interval as minijets. For dijet 
events at large rapidity intervals the BFKL theory systematically resums the leading 
powers in the rapidity interval y, including both real and virtual gluon corrections. 

In ref. |^ we showed that the exponential enhancement with y in dijet production at 
fixed xa and xb , originally suggested as a signature of the BFKL minijets by Mueller 
and Navelet0, is highly suppressed by the parton distribution functions at Tevatron 
energies. However, we found that other observables such as the transverse momentum 
p± distribution and the jet-jet correlations in p± and azimuthal angle are significantly 
affected by the minijet resummation. For example, we saw that these correlations are 
not a leading feature of the expansion in the rapidity interval. Accordingly, they fade 
away as the rapidity interval increases. 

However, in comparing the BFKL resummation with the exact 0{al) calculation^], 
we noticed that the large-rapidity approximation to the kinematics causes a serious error 
in the BFKL predictions when the tagging jets are not back-to-back in p± and 0. In 
order to account for this error we introduced an effective rapidity y which restricts the 
phase space of the minijets in such a way that the truncation of the BFKL resummation 
to 0{a^) agrees with the exact 2^3 0{a^) calculation. For large y, the difference 
y — y is nonleading. Since the rapidity variable which is resummed by BFKL is only 
defined up to transformations y ^ y + X where X is subleading at large rapidities, we 
used y instead of y in the BFKL resummation in order to obtain quantitatively more 
reliable predictions of the transverse momentum distributions. We use here the effective 
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rapidity y in the BFKL resummation to analyze the moments [|T0| of the decorrelation in 
azimuthal angle. 



2 The minijet resummation and the effective rapid- 
ity interval 

We are interested in the semi-inclusive production of two jets in hadron-hadron collisions 
PaPb jij2+X. We describe the two partonic tagging jets by their transverse momenta 
and rapidities and (p2±,l/2), where we always take yi > y2- We reexpress 

the jet rapidities in terms of the rapidity interval y = yi — y2 and the rapidity boost 
y = {Vi + l/2)/2. This is convenient since we are mainly interested in the behavior of the 
parton subprocess, which does not depend on y. For large values of y the cross section 
for this process can be written 



dy dy dpi^dp2±_d4> dpij_dp2_i_Lii^ 

where the parton momentum fractions are dominated by the contribution from the two 
tagging jets 

_ Pi±e^' 

4 ^ (2) 

and fi is the factorization/renormalization scale. In this limit the amplitude is dominated 
by gg, qg, and qq scattering diagrams with gluon-exchange in the t-channel. The relative 
magnitude of the different subprocesses is fixed by the color strength of the respective jet- 
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production vertices, so it suffices to consider only gg scattering and to include the other 
subprocesses by means of the effective parton distribution function fes{x,fi'^) [lTl|,|f^. 

The higher-order corrections to the gg subprocess cross section in (|I|) can be expressed 
via the solution of the BFKL equation0, which is an all-order resummation in of the 
leading powers of the rapidity interval 



(3) 



with 



and 



M4i/,pi^,P2±) = rciz/e-("''^)^cosfz/ln4^') 



(4) 



u;{n, v) = ^^[^(1) - Re^(^^ + (5) 

TT Z 

and ip the logarithmic derivative of the Gamma function. Eq. (j^) can be expanded order 
by order in as- At 0{as) and for pi± ^ p2± we obtain 



dpl^dpl^dcf) ^t^pI^ pI^ pI^ + pI^ + 2pi_LP2± cos ' 
This can be compared with a fixed-order calculation of dijet production at the same 
order of computed through the 2^3 parton amplitudes, 



dydydplj_dplj_d(j) Iy2 ^^'^ I ^^^^^^''"^^^^^ ^ ^ ^ dpl^dpl^dy^dcj) 



where is the rapidity of the third final-state parton, integrated over the interval 
spanned by the tagging jets. = Q,Q,G labels the distribution function of the 

parton species and flavor = q,q,g inside hadron A{B). We include all parton 
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subprocesses ||12|| , and use the exact values of the parton momentum fractions 



xa = 



(8) 



Xb = 



In order to avoid the colhnear singularity between final-state partons, which is however 
subleading at large rapidities||], configurations in (0) where the distance R between two 
of the partons on the Lego plot in azimuthal angle and rapidity is smaller than the jet 
cone size Rcut are discarded. In the large-y hmit (|^) reduces to (|1|), with the parton cross 
section given by (||). 

However, we have seen in ref. [|^ that the large- ?/ approximation seriously overesti- 
mates the cross section when the two tagging jets are not back-to-back in pj_ and 0, even 
for rapidity intervals as large as ?/ = 6. This occurs because the large-y cross section (^) 
assumes that the third (minijet) parton can be produced anywhere within the rapidity 
interval [y2,yi] with equal probability, whereas in the full 2 — 3 cross section the prob- 
ability is highly suppressed by the structure functions when the third jet strays too far 
from the center of this interval. In ref. we introduced an "effective" rapidity interval 
y to take into account the fact that the range in rapidity spanned by the minijets is 
typically less than the kinematic rapidity interval y. We define y{n,pi±,p2j_,y,y) by 



where n is the Fourier series index of eq. (|]). The cross section in the numerator is that 
of eq. (1^) and is computed using the exact kinematics (1), while the cross section in the 
denominator is that of eq. (|^) and is computed using the large- ?/ kinematics (0). The 
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denominator can easily be computed analytically using the large-y solution (||). Note 
that y is defined so that if we replace y —>■ y in the BFKL solution @ and truncate to 
0{ag) we recover the exact 2^3 cross section. 



3 The azimuthal angle decorrelation 

Our aim is now to reconsider the azimuthal angle decorrelation H],|TO[ by using the 
effective rapidity y instead of the kinematic one y in the BFKL resummation. Using y, 
we expect that the distribution will be more populated around the peak (p = tt where 
the tagging jets are back-to-back, and less populated on its tails where the tagging jets 
are decorrelated, as compared to the decorrelation plot we considered in ref. [§]. This is 
because the larger the decorrelation, the larger the overestimating error in approximating 
the parton momentum fractions with eq. (|^). However, in practice the computation of 
the Fourier coefficients @ using eq. (H) for y becomes quickly very time-consuming as 
n grows. To test the BFKL picture, it is then more convenient to compute the moments 
of the azimuthal angle [O], defined by 



Jq"" d4) cos n{4)-7i){dao/ dyd(j)) 

< cosnlcp — n) > = — . (10) 

/q d(j){dao/dyd(j)) 

For a 5-function distribution at = tt, as occurs at the Born level, all of the moments 
will equal one, while for a fiat distribution all of the moments will equal zero for n > 1. 
Thus, the decay of the moments from unity is a good measure of the decorrelation in (p. 

In Fig. 1 we plot in the solid curves the first two moments of the azimuthal angle 
as a function of y, calculated using the effective rapidity y in the BFKL solution. We 
use the LO CTEQ2 parton distribution functions ||13|| with the ren./fact. scales set to 



/^^ = pi±p2±- To facilitate the comparison with the experimental data we choose the 
kinematic parameters of the data analysis of the 1993 run at the DO detector. Namely, 
we require that one of the tagging jets must have p± > 50 GeV, while the other must 
have p± > 20 GeV. The rapidity boost y is integrated over, subject to the constraint 
|z/i|maa; = \y2\max = 3.2. The rapidity interval is integrated in unit bins centered around 
the variable y of the plot. Finally, to compare with our previous analysis, we also show 
in the dashed curves the moments calculated using the kinematic rapidity. 



1.0 




Figure 1: The moments of the azimuthal angle distribution for n = 1 and 2, computed 
with the BFKL resummation using the effective rapidity y (solid) and using the kinematic 
rapidity y (dashed). The error bars are the estimated statistical error from the monte 
carlo integration. 



As expected, the naive BFKL analysis (dashed) tends to overestimate the decorrela- 
tion in azimuthal angle. The first two moments are typically two to three times smaller 
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than those calculated using y in the BFKL solution (solid). However, there is still a 
substantial decorrelation exhibited in the solid curves which should be experimentally 
measurable. For small rapidities one would expect these results to be invalidated by 
nonleading effects, but for larger rapidities (optimistically, perhaps y ^ 4) the minijet 
effects should begin to dominate. This is a particularly nice process for probing the 
BFKL structure, from both an experimental and a theoretical point of view. Experi- 
mentally, the measurements are much less affected by detector resolution than the p_L 
measurements . Theoretically, the azimuthal angle distribution is less sensitive to the 
factorization/renormalization scale used. Thus, a measurement of the azimuthal angle 
decorrelation should be an excellent probe for the validity range of the BFKL dynamics. 
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